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Inflation in a five-dimensional brane world model with two boundary branes is studied. We 
make use of the moduli space approximation whereby the low energy theory reduces to a four- 
dimensional biscalar-tensor gravity plus a minimally coupled scalar field. After a detailed analysis 
of the inflationary solutions, we derive the evolution equations of the linear perturbations separating 
the adiabatic mode from two entropy modes. We then examine the primordial scalar and tensor 
power spectra and show that their tilt depends on the scalar-tensor coupling constant. Finally, the 
induced CMB anisotropies are computed and we present a Monte Carlo Markov Chains exploration 
of the parameter space using the first year WMAP data. We find a marginalized probability bound 
for the associated Eddington parameter at the end of inflation 1 — 7 < 2 x 10~ 3 , at 95% confidence 
level. This suggests that future CMB data could provide crucial information helping to distinguish 
scalar-tensor and standard inflationary scenarios. 

PACS numbers: 04.50.+h, ll.10.Kk, 98.80.Cq 



I. INTRODUCTION 

Although the standard particle physics and cosmo- 
logical models in a four-dimensional spacetime provide 
an accurate enough description of particle physics ex- 
periments and cosmological observations p|, consider- 
ing unseen extra dimensions allows one to make theo- 
retical progress towards a Unified Theory of all interac- 
tions El El El IS 0- The advent of M-theory and branes 
has led to the "brane world" concept in cosmology. Moti- 
vated by the Hofava-Witten model and its compactifica- 
tions |a,|9|i brane world models assume that we are living 
on a higher dimensional spacetime boundary (brane). 

In cosmology, brane worlds have been studied within 
the framework of five-dimensional (5D) models and lead 

to only a few observable predictions (see Refs mm [2 

for a review and references therein). Models ranging 
from classical membrane theories of General Relativity 
to topological defect models have been studied in vari- 
ous compact or non-compact, curved or flat, symmetric 
or non-symmetric, extra dimensional spacetimes, as well 
as their associated deviations, or incompatibilities, with 
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res pect to the standard cosmolo gy lie 

p3 p| p| EE EI iSha i2fj, im 

I34L l3fl ISfl l37l l88| . Definite cosmological predictions re- 
main however difficult to extract in the high energy limit 
where the universe evolution is strongly affected by the 
physical processes occurring in the extra dimensions. On 
the other hand, at lower energy scales (typically smaller 
than the brane tension) it appears possible to describe 
such systems by effective four-dimensional actions. 

In the framework of five-dimensional compactifica- 
tions of M-theory [39L l4fj| . the moduli space approxi- 
mation describes, through a four-dimensional effective 
action, a system of two branes of opposite tension em- 
bedded in a five-dimensional warped spacetime |4l| . 
Besides the fields living on the positive tension brane 
(assumed to be our universe), the moduli associated 
with the position of the branes in the fifth dimension 
act as two non-minimally coupled scalar fields thereby 
leading to an effective biscalar-tensor theory of grav- 
ity lIllIllIllIllIilllE Scalar-tensor theo- 
ries have been intensively studied as natural extensions 
of General Relativity and are strongly constrained in the 
solar system [5(3, |^. On the cosmological side, there 
are (weaker) constraints coming from the Big Bang Nu- 
cleosynthesis (BBN) and recent works have been focused 
on scalar-tensor effects in the Cosmic Microwave Back- 
ground anisotropies (CMB) and in weak gravitational 
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these works, scalar-tensor theories have been discussed in 
the context of the post-inflationary cosmology by quanti- 
fying the distortions the non-minimally scalar fields may 
induce compared to the pure General Relativity predic- 
tions. In fact, even if the cosmological constraints are 
currently weaker than the solar-system ones, they are 
complementary in the sense that they probe the cosmic 
times. This is relevant since it has been shown that 
scalar-tensor gravity models can relax towards General 
Relativity during the cosmological evolution [(Ml. As a 
result, such models may have no observable impact in 
the post-inflationary eras, and in particular in the solar 
system, while modifying significantly the dynamics of the 
early universe. 

On top of probing the post-inflationary eras, the preci- 
sion reached by the current CMB experiments allows one 
to constrain the shape of the power spectra of the primor- 
dial cosmological perturbations. Recast in the inflation- 
ary context, this is an opportunity to shed light on the 
physics at earlier times and in particular the deviations 
scalar-tensor theories may generate during inflation. 

In this paper, we are interested in the primordial in- 
flationary eras occurring in a boundary inflation model 
and their potential observable effects in the CMB, given 
the first year Wilkinson Microwave Anisotropics Probe 
(WMAP) data [HI El E3- A minimal realistic three- 
fields model is considered: two of them being the mod- 
uli associated with the position of the branes in the ex- 
tra dimension, and the other acts as a standard infla- 
ton field li ving on our brane and coupled to the matter 
sector 144L |45L liti, l6Sj . After a detailed analysis of 
the several viable background evolutions, we compute 
and discuss the primordial power spectra of the scalar 
and tensor perturbations generated during inflation. We 
then derive the resulting CMB anisotropics and perform 
a Monte Carlo Markov Chain exploration of the parame- 
ter space in a regime where the deviations from General 
Relativity predictions are dominated by the shape of the 
primordial power spectra. By comparing the power spec- 
tra induced by the standard chaotic inflaton field alone 
and the ones obtained by considering the moduli, we find 
that the first year WMAP data can differentiate between 
standard chaotic and boundary inflation. In fact, for a 
given matter sector, the presence of non-minimally cou- 
pled scalar-fields during inflation changes the tilt of the 
primordial power spectra due to the running of the con- 
formal factor. In our framework, we find the presence of 
the moduli statistically disfavored by the data, either as 
the WMAP data are insensitive to their effects when they 
are weakly coupled, or as they lead to a strongly tilted 
power spectra. These effects lead to an upper bound on 
the allowed values of the moduli coupling constant during 
inflation, which can be recast in terms of the more com- 
monly used Eddington post-Newtonian parameter |5l| : 
1 — 7 < 2 x ICU 3 . On one hand, this bound holds at the 
time of inflation and, according to the above discussion, 
provides a very early constraint on scalar-tensor theories. 
On the other hand, it relies on the running of the con- 



formal factor during inflation as well as the choice of the 
matter sector. In particular, one may expect this bound 
to be modified, or evaded, with more complicated poten- 
tials for the matter fields, or by freezing the running of 
the conformal factor at the time where observable cosmo- 
logical perturbations were generated. Nevertheless, this 
result suggests that the current and future CMB data 
can be used to probe scalar-tensor theories as early as 
during the inflationary era. 

In Sect.[n] the boundary inflation model is introduced 
and the background inflationary evolutions solved both 
analytically and numerically. Sect. IIIll is devoted to the 
derivation of the linear perturbations around the back- 
ground evolution by means of a separation between the 
adiabatic and the two entropy modes. The equations of 
motion obtained arc solved numerically and used to com- 
pute the resulting primordial power spectra at the end of 
inflation. In Sect. IIVI we use modified versions of the 
public CMB codes CAMB |g<| and COSMOMC linked 
with our inflationary code, to probe the boundary model 
parameter space given the WMAP data. Finally, the rel- 
evance and the possible improvements of our results are 
discussed in the conclusion, and the possibility of using 
the CMB to accurately probe scalar-tensor theories in 
the early universe is raised. 



II. BACKGROUND EVOLUTION 

We will focus on the low energy effective theory of 5D 
brane world models with a bulk scalar field and matter 
living on the positive tension brane. In this context, the 
low energy degrees of freedom of the model are the gravi- 
ton leading to 4D gravity, the bulk scalar field zero mode 
and the inter brane distance, i.e. the radion. One can 
also parametrise the two scalar modes as the two brane 
positions. In the following we will denote by ip and tfj the 
two moduli which are related to the brane positions z± 
and z 2 by 



c v coshi/; = (l - 4fc b c^z 1 
e v sinh ip = (l — 4fcbC 2 v z 2 



(2 C 2+l)/4c" 



(1) 



where fcb fixes the energy scale of the brane tensions and 
c w is a constant which determines how warped the ex- 
tra dimension is |4fLl4Ml4r|. In particular for c w = 0, we 
obtain the Randall-Sundrum case. The four-dimensional 
(4D) gravitational coupling is related to the 5D gravita- 
tional coupling via 



K 2 = (24 + I) k h K 2 5 . 

Using these relations one gets 



[1 - 4fc b c w z 1 ) 



(l-4fc bC ^ 2 ) (2c " +1)/2cS 



(2) 



(3) 
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and 



tanh^ — 



1 - 4k h cjz 2 
1 - 4fc b c 2 v z 1 



(24 + l)/4<4 



(4) 



Notice that <p — > — oo and ip — > oo when the two branes 
collide. On the contrary, ip — > and y> — ► oo corre- 
spond to the negative tension brane being stuck at the 
bulk singularity z 2 — and the positive tension 

brane receding away towards infinity. This situation will 
happen during the inflationary evolution. 

Denoting by \ the minimally coupled scalar field living 
on our brane universe, the effective action in the Einstein 
frame, with a metric of (— , +, +, +) signature, reads 



S = 



1 

2k 5 



R - c p (dip) 2 - c r (dip) 2 - W V=gd 4 x 



A 4 U 



-gd A x, 



(5) 

where tp and ip stands for the moduli fields, A(ip, ip) is 
the conformal factor, W((p, ip) the moduli bulk potential 
and U (x) the minimally coupled scalar field potential on 
the brane. In the following, we will focus on the minimal 
setup of Refs. |4jj, i.e. 

A 2 (ip, ip) = e ~( c p/ 3 ^ (cosh iP) Cl/3 , W((p, ip) = 0, 

12c 2 , 6 (6) 

Cr 



2k 2 



R - W" - 2V(T C ) 



-gd x, 
(9) 



the field-manifold being defined through 

T a = [ i I , Lb =diag(^ 2 ) c p ,c r ) , (10) 

where the dimensionless scalar field x — K X has been in- 
troduced. From Eq. iJBJ, the fields evolve in the potential 



V = k 2 A 4 U. 



(11) 



In the following, it will be implicitly assumed that the 
Latin indices {a, 6, c, d} refer to the field-manifold. Differ- 
entiating the action @ with respect to the metric leads 
to the Einstein- Jordan equations 



G in, S in, , 



with the source terms 



(12) 
(13) 



where 



Sf v = d^ a d v T h - -g^d p ^ a d^ b . (14) 



Similarly, the fields obey the Klein-Gordon-like equation 



In terms of the 5D picture, c w parametrises the coupling 
of the bulk scalar to the brane. Since the resulting four- 
dimensional theory is of multiscalar-tensor kind, it is also 
convenient to introduce the first conformal gradient 



din A 

dip 



dlnA Cr 

^/- = -tanhV>. (7) 



These are the gravitational couplings of the two mod- 
uli. Their present values are constrained by solar system 
experiments. Moreover, we will assume a chaotic-like po- 
tential for x 



U{x) = \m 2 c x 2 . 



(8) 



In such a setup, the field x 1S expected to be coupled to 
the observed form of matter and it will be referred to as 
the "matter field" in the following. As in the standard 
chaotic model, inflation can be driven by x, but also by 
the moduli as discussed in the next section |tT| . 



A. Basic equations 

1. Sigma-model formalism 

For the sake of clarity it is more convenient to recast 
the action |J5J in terms of a non-linear sigma-model [Hcl 



UT c + g^^g^ag^b = yc ^ 



(15) 



where T denotes the Christoffel symbol on the field- 
manifold 

Tafc = ^ Hda.b + U,a ~ Lb,d) , (16) 

and V c should be understood as the vector-like partial 
derivative of the potential 



yc = £ cdy d _ £ cd 



dV 

dlF 3 ' 



(17) 



2. Equations of motion 

In a flat Friedman-Lemaitre-Robertson- Walker 
(FLRW) Universe with metric 

ds 2 = g^dx^dx" = a 2 { v ) {-drf + * y da?W) , (18) 

t] being the conformal time, the equations of motion l|12|) 
and l|15fl simplify to 



m 2 = -£abF a 'F b ' + a 2 V, 



2H' + K z = -^ab^ 1 '^' + a 2 V, 



(19) 
(20) 



T c " + T c ah T a 'T hl + 2UT cl = -a 2 V c , (21) 
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where a prime denotes differentiation with respect to the 
conformal time and TL is the conformal Hubble parameter 
H = a! I a. In the Einstein Frame, the "new terms" com- 
pared to minimally coupled scalar multifield inflationary 
models are encoded in the sigma-model metric £ a b and 
the Christoffcl symbols T^ h . 

Since we are interested in inflationary behaviour, it 
is more convenient to work in terms of an "efold" time 
variable 



hi 



(22) 



Apart from being the relevant physical time quantity dur- 
ing an inflationary era, when expressed in terms of n, the 
field equations can be decoupled from the metric evolu- 
tion Hfl [zl III El and Eqs. JTgJ) to (JUJ separate 
into 



H 

H 
H 



V 



3- ;<r 2 



T c = - — 



XI 
v 



(23) 

(24) 
(25) 



The dot denotes differentiation with respect to n, the 
physical Hubble parameter is H = H/a, and the field 
velocity (in efold time) on the manifold is 



(26) 



This is the derivative of the so-called adiabatic field. In- 
deed, differentiating Eq. Q24JI and using Eq. 1)2 5J1 yields 
the standard second order differential equation for the 
adiabatic field <r, i.e. in conformal time [76L l77| 



cr" + 2Ha' + a 2 L c V c = 0, 



(27) 



where t a are unit field-vectors along the classical trajec- 
tory 



while the only non-vanishing sigma-model Christoffel 
symbols are 



xv 



xx 

XX 



-A 2 a^, 



-A 



(30) 



where lowering and raising indices on the first conformal 
gradients is performed with the field- manifold metric £ a b. 
The fact that the effective potential for the fields is ln(V^) 
ensures that the evolution of the fields in efold time is in- 
dependent of any multiplicative constant appearing in 
the definition of the matter potential U 7 as m c in the 
chaotic- like potential we are interested in [see Eq. ©]. 
Obviously, m c still fixes the normalisation of the Hub- 
ble parameter in Eq. I|23[) and is directly related to the 
amplitude of the power spectrum of linear perturbations. 

From Eqs. (|25|l and (|29|l . we can draw the qualitative 
evolution of the fields. Starting from reasonable initial 
conditions, i.e. a 2 < 6 [see Eq. i)23[) ]. the friction terms 
will tend to produce a stationary regime in which the 
fields relax toward the minimum of the potential ln(l^). 
For instance, the dilaton field ip will relax toward vanish- 
ing value as fast as 



—n. 

3 ' 



(31) 



and stays frozen at -0 = afterward, whereas Lp grows as 



(32) 



From Eq. (|29|) . the matter field \ is expected to slowly 
evolve from its initial value provided 4 2 y> 2, thereby 
allowing a period of slow-roll inflation [43|. However, 
since the force term involves also the conformal factor 
A 2 , the properties of the inflationary era, as its duration 
and end, will be certainly dependent on the evolution of 
the moduli fields tp and ip. 

In the next section, we will present some analytical 
approximations and numerical solutions to the previous 
sketch. 



T" 



& 



As mentioned above, the fields evolution is decoupled 
from the Hubble flow and Eq. (|25|l mimics the evolu- 
tion of a "particle" of variable inertia 1/(3 — fr 2 /2) on a 
curved manifold in presence of friction and an external 
force -V c /V Hfllfil- From Eqs. © and © one gets 



ye 2 
— = 4a v = -- 

V 3' 

i 2 

= 4a v = - tanhi/>, 

V 3 v ' 



i u* 



1 2 



(28) B. Inflationary solutions 

1. Scale factors acceleration 

Since we are dealing with a multiscalar-tensor theory, 
it is convenient to clarify the definition of "inflation" in 
the cosmological context. The frame where non-gravity 
measurements take their usual form is the Jordan frame 
where matter is minimally coupled to the metric. For a 
given conformal coordinate system, the metric tensor g 

(29) m th° Jordan frame is obtain from the metric tensor g in 
the Einstein frame by the conformal transformation 



A 2 U A 2 x 



9 = A 2 g 



Aa. 



(33) 
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As a result, the acceleration of the scale factor is not 
necessarily the same in the Jordan and Einstein frame 
and may even take place in one frame only jj^ • To solve 
the flatness and homogeneity problem, we are interested 
in inflationary era in the Jordan frame, i.e. when the first 
slow-roll parameter in that frame 



f 



1<i. 

u 2 



(34) 



From Eq. I|33[l . the conformal Hubble parameter in the 
Jordan frame is also 



H = H(l + a a T a , 



(35) 

where the dot denotes again the derivative with respect 
to n, the total number of efold in the Einstein frame. 
From Eq. I|34f> . one gets 

. e + a a T a a a T a + fUF a F b 



1 + a a T a 



1 + a a T a 



(36) 



where e is the first slow-roll parameter in the Einstein 
frame and /3 a b is the second conformal gradient defined 
as 

da a 



flab 



(37) 



In our setup, from Eqs. © and 0, its only non- vanishing 
value is 

Aw = ~K TT7' ( 38 ) 

o cosh ip 

which rapidly vanishes for non- vanishing values of ip. 
Moreover, in a friction dominated regime T a ~ 0, the 
remaining term in Eq. (|36|l can be significant only at the 
time when both (3^ and ip are non-vanishing, i.e. when 
the ip field leaves the friction dominated regime, where it 
was decreasing with constant velocity, to be frozen on the 
attractor at ip = 0. As the result, we can expect e ~ e 
as long as the fields remain in the friction dominated 
regimes. 



2. Minimal slow-roll approximations 

From the Einstein equation l|24[l . the first slow-roll pa- 
rameter in the Einstein frame is simply related to the 
field velocity by 

e = = \A*j? + + \c T ip\ (39) 

Putting everything together in Eq. l]25f) yields 

/ c c \ . 2 

X = - (3- e- yV? + yV'tanh-0J X - ^rr(3 - e), 



A 2 x 



Cp = -{Z-e)v-\A 2 X 2 + \{Z-e), 



(40) 
(41) 



i> = - (3 - e) ip + tanh^) - -(3 - e) tanh^. (42) 
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The friction dominated evolution of the moduli fields 
discussed in the previous section is recovered provided 
we assume the matter field to be in slow- roll Ax <C 1. In 
this limit, Eqs. (|41|> and (|42[) admit the solutions 



ip 



2 
3' 



(43) 



where it was assumed that ip > 1. However, note that 
-0 = is an exact solution of Eq. I|42|) and thus, ip may 
be non-vanishing only during less than a few efolds of 
transition between ip ~ 1 and ip — 0. Let us point out 
that the derivation of Eq. (|43|l docs not require all the 
fields to be in slow-roll, and in particular e may not be 
small [?3 ■ This is due to the simple form of the conformal 
factor we are interested in, whereas in the gen eric case 
a full slow-roll treatment may be required |80L l8l| . In 
the two regimes ip > 1 and "0 = 0, the conformal factor 
evolves as 



A 2 



A 2 p- Cn 



(44) 



with C = 2(c r + c p )/9 for ip > 1, and C = 2c p /9 for 
ip = 0. In both cases, since we are assuming Ax <C 1, we 
have e ~ C [see Eq. I|39|l ]. The moduli are thus driving 
two expansion eras with 



e - g (cr + c p ) (V > 1), 



-c p (ip = 0). (45) 



In our setup, according to Eq. ©, the first era is not 
inflation since e ~ 4/3 > 1, whereas the second one is 
for values of the coupling constant c p < 9/2. Physically, 
it means that starting from high values of ip only gives 
rise to a period of non-inflationary expansion where the 
moduli ip relaxes toward ip = 0, and consequently, the 
effective Planck length K e g — An decreases exponentially. 
When ip reaches its attractor, an inflationary era driven 
by ip can start, which also corresponds to a decreasing 
effective Planck length [see Eq. (|4*4l ]. 

In terms of the 5D picture, the fate of the moduli in- 
dicates that the two branes move away from each other. 
The inflationary phase driven by ip occurs once the right 
(negative tension) brane becomes close enough to the 
bulk singularity. As x i s constant in this regime, the 
energy density of the boundary inflaton is constant lead- 
ing to a constant detuning of the left brane tension. As 
already studied in [4l| , this leads to an accelerated phase 
with an exponential potential and power law inflation. 
More details about this power law inflation phase can be 
found in Ref. plf. 

This picture remains valid provided Ax <C 1 and from 
Eqs. (|40|l and (|43|) we can check the consistency of this 
approximation. From Eq. (|44H . if x starts from large 
enough values such that AiaiXm 3> \/2/ e then an ap- 
proximate solution of Eq. (|40|l is 



X~Xh 



eA Lxv 



(46) 
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In this case, we have indeed 



^iniXi] 



«1, 



(47) 



and the matter field \ remains frozen at its initial value 
until eA 2 \i n i — 1, i.e. during approximately 



1 



An ~ - In (A 2 xf ni e) 



(48) 



cfolds after which Eq. (|46[) is no longer valid. Then the 
field x starts to evolve and A\ increases until e = 1 for 
which inflation ends. Nevertheless, the standard slow-roll 
approximation can be used to derive the behavior of the 
matter field during this last phase. For reasonable small 
value of c p , e ends up being dominated by the matter 
field evolution [see Eq. whereas the behavior of the 
moduli remains almost the same. Indeed, as can be seen 
in Eq. I|41(l . the value of <p will be affected up to 10% only 
by the matter field when A\ — v2, which is also the end 
of inflation [see Eq. l|39"j) ]. As a result, Eq. lf4*4"jl is still a 
good approximation for the conformal factor, and using 
the slow-roll approximation x <C x in Eq. (|40|l leads to 



A 2 



C 



^Cn sr 



c 



where 

1 2 -62 



the "sr" indices label the efold at 



A 2 x 2 /2. The end of inflation occurs for A 2 x 2 
gives the number of efolds in this last stage 



(49) 

which e ~ 
2 which 



C 



In 



V 



- —A 2 ■ y 2 

^ miAsr 

c 

1+ 2 



(50) 



Note that in the limit C — > 0, we recover the usual one- 
field slow-roll expression for a chaotic potential js3 



A; 



1 2 2 
7 K cffX sr 



(51) 



with K c ff = A- Y 



3. Numerical approach 



A numerical integration of the equations of motion ll'MI) 
to 125fl has been performed and illustrates the analytical 
behavior described in the previous section. In Fig. ^ 
the three expansion eras can be distinguished by means 
of the e values. For the chosen model parameters, dur- 
ing the first ten efolds, the field if) relaxes towards zero 
and its evolution dominates the expansion. As previously 
pointed out, this is not an inflationary era in our model. 
During this phase, the conformal factor, and thus the 
Hubble parameter, decrease exponentially [see Eqs. i|ll|) 
and 1(23)) ]. 

Once if) reaches its minimum, the dynamics are dom- 
inated by the second moduli ip and e ~ 2c p /9. Note 
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FIG. 1: Typical behaviors of the matter field \ an d the two 
moduli ip and if) plotted with respect to the total number of 
efold. The Hubble parameter in unit of the matter field bare 
mass m c , and the conformal factor, as well as the first slow- 
roll parameter in Einstein and Jordan frame are shown. As 
seen from the e behavior, this is a three stages expansion. The 
model parameters are c p = 0.01, c r = 6 — c p [see Eq. ||SJ], while 
the initial conditions have been chosen such that ip\ n i — 4.5, 
Ai„i = 1000 and Xini = 5. Moreover, the field velocities have 
been chosen to start immediately in the friction dominated 
regime, i.e. with T a = -V a /V (see also Fig.^J. 



that during the transition between these two regimes, 
the slow-roll parameter e and e, in the Einstein and Jor- 
dan frames respectively, do not match, as expected from 
Eq. J2HJ). As can be seen on Fig. ^ inflation (in the 
Jordan frame) starts a few efolds before e < 1 in the Ein- 
stein frame. During the p> dominated inflationary era, 
the matter field x leaves its initial value and starts to 
slow-roll under the external force — V*/V [see Eq. I|29ll ] 
leading to the observed variation of e after n ~ 100. The 
last stage occurs when the inflation is mainly driven by 
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efold 

FIG. 2: Evolution of the fields for a two stages expansion. The 
matter field \ is destabilized as soon as ij) reaches its attractor 
tp = 0. As the evolution of e shows, there is not time for a 
pure tp dominated inflationary era and we start immediately 
in a x slow-roll dominated inflationary era until e = 1. These 
solutions have been obtained for the same model parameters 
as in Fig. Q i.e. c p = O.Of, and for the initial conditions 
tp ini = 4.5, A ini = 400 and Xmi = 2. 



X and we recover a one-field like behavior until e = 1 for 
which inflation stops. Note that, as expected from the 
analytical study, even in this last regime, the evolution 
of the moduli tp and ip is not significantly affected by 
the rolling of the matter field. As a result, their second 
derivatives remain small and e ~ e at the end of inflation 
(see Fig. HJ. The above settings can nevertheless be 
affected according to the parameter values and the ini- 
tial conditions. Indeed, the matter field \ can leave its 
initial value as soon as, or even before, the end of the tp 
domination expansion. In that case, as can be seen in 
Fig. |2 there is no time for a pure tp dominated inflation 
and after a intermediate mixed phase, we jump directly 
to the last x dominated stage. 

Concerning the initial conditions, it is important to 
mention that we have fixed the initial derivatives of 
the fields in such a way that their evolution start in 
the friction dominated regime, i.e. with an initial boost 
jra = _v a /V. In Fig. El we have plotted the behavior of 
the fields obtained for an arbitrary choice of reasonable 
initial velocities, i.e. verifying a 2 < 6. As expected, it 
takes less than few efolds for the fields to relax toward 
their attractor behavior. As previously discussed, during 
the relaxation time, the slow-roll parameters in Einstein 
and Jordan frame are significantly different. 

In order to derive the primordial perturbations pro- 
duced in this model and the resulting cosmic microwave 
background anisotropies (CMB), it is necessary to as- 
sume a cosmological scenario for the background. Mo- 
tivated by the brane world picture where the moduli tp 




„ 

0.01 0.1 1 10 100 

efold 



FIG. 3: Relaxation toward the friction dominated regime. 
The dash and thick curves represent the evolution of the field 
velocities (in efold time) and the slow-roll parameter for initial 
conditions such that tp- m i = 0, Aj n i = 1, Xmi = 16, c p = 0.01, 
with arbitrary chosen initial field velocities: xim = 1, <Pini = 
—0.5 and i/'im = 0.6. The dotted curves are obtained for 
the same model parameters and initial field values, but with 
initial velocities on the attractor T a = — V a /V. As expected, 
it takes less than a few efolds for the velocities to relax toward 
their friction dominated values. Note the difference between 
e and e during the relaxation time. 
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e e (Einstein frame) 










e e (Jordan frame) 




















i , 
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efold 

FIG. 4: The end of inflation and the beginning of the re- 
heating phase where the matter field oscillates. Note the new 
slope tp ~ 1/3 coming from the average value of the slow-roll 
parameter during that stage [see Eq. I|41|l ]. Since tp and tp are 
not dramatically affected by the oscillations, one still have 
e ~ e. 
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and ip are related to the position of the branes in the 
five-dimensional spacetime, it is reasonable to consider 
that these two fields do not decay and remain present 
in the late-time cosmology. As a result, only x can de- 
cay into radiation and dark matter at the end of infla- 
tion. In Fig. wc have plotted the fields evolution dur- 
ing a few efolds after the end of inflation, the matter 
field x ends up oscillating thereby starting a period of 
reheating from which the radiation dominated era can 
begin [83l l84l l85j . Note that even if inflation can be 
driven by ip (see Sect. Ill B|) . the presence of x is still re- 
quired to end it. The observational consequences coming 
from these assumptions will be more detailed in Sect. lIVI 
In the next section, the linear perturbations arising 
during the inflationary eras are discussed, as well as their 
primordial power spectra. 



III. LINEAR PERTURBATIONS 



Gravitational and matter sector 



The fact that there is more than one scalar field in- 
volved leads to the existence of entropy modes that can 
source the adiabatic mode. In the next section, we use 
the formalism developed in Refs. [7(| 173, EE 113 to de- 
rive the evolution equations of the adiabatic and the two 
entropy modes arising in our model. 



B. Adiabatic and entropy perturbations 

1. Generic decomposition 

From Eqs. to (JSSJ, up to the background equa- 
tions, the evolution of the Bardeen potential simplifies to 



+ eft*' + (2H' + AH 2 ) * - A* = -a 2 V c 5T c . (58) 



In terms of the comoving curvature perturbation 88, 8i 



In the longitudinal gauge, the scalar perturbations 
(with respect to the rotations of the three-dimensional 
space) of the FLRW metric can be expressed as 

ds 2 = a 2 [- (1 + 2$) drj 2 + (1 - 2f) 7 ^dxMa; J '] , (52) 

where the indices i and j refer to the spatial coordinates 
only, and <& and '5 are the Bardeen potentials. In the 
sigma-model formalism, denoting by bT the scalar per- 
turbations of the fields, the Einstein- Jordan equations 
perturbed at first order are 



3W + (H' + 2H 2 ) * - A* = -^tabT a '8T hl 
- ~ (-Lb^'F 1 " + a 2 v)j 5T C 



#" + 3H^' + {W + m 2 ) W = t zi a hT a 'bT b ' 



1 (1 



2 V 2 



l*^cT a 'T h ' ~ a 2 V c ) 6F 



(53) 
(54) 

(55) 



where the prime stands for the derivative with respect 
to the conformal time and use has been made of $ = 
(perturbed Einstein equation for i ^ j). The Laplacian 
reduces to 



A = S lJ didj, 



(56) 



for flat spatial hypersurfaces. The perturbed equations 
of motion for the fields stem from Eq. I|15|) and read 

5T C " + 2T L ab T al 5T hl + 2H5T cl 

- AST C = m'T cl - 2^a 2 V c . (57) 



n 



w -n 2 



Eq. JI2J yields 



C = * + H% 



(59) 



(60) 



where the adiabatic perturbation da is the perturbed ver- 
sion of Eq. I|26|l as well as the resulting perturbation 
of all fields projected onto the classical trajectory [see 
Eq. JUJ]: 



Zab-r OJ- 
0(7 = ; = L a bT . 



The dynamical equation (|58l) also reads 



(61) 



2H 



2H 



C = ^A* - ^ ( a 2 V a 8T a - a 2 



which can be recast into 



V C T CI £ ab J ra, 5J rb 



C = ^A* a 2 V c 5T d . 



(62) 
(63) 



O 0" 

The orthogonal projector is defined by 

-Ub = Lb ~ Vab, (64) 

where rj a b = ta^b is the first fundamental form of the 
one-dimensional manifold defined by the classical tra- 
jectory [94|]. We recover that the comoving curvature 
perturbation on super-Hubblc scales (A\I/ ~ 0) is only 
sourced by the entropy perturbations defined as the pro- 
jections of all field perturbations on the field-manifold 
subspace orthogonal to the classical trajectory [9^ . 
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2. Spherical basis 

For the studied boundary inflation model, the sigma- 
model manifold is three-dimensional which implies the 
existence of two entropy modes living in the two- 
dimensional subspace orthogonal to the classical trajec- 
tory. The decomposition performed in Refs. [7(1 I77| is 
straightforwardly generalized by choosing a local spheri- 
cal basis at each point of the fields trajectory. This can 
be performed by introducing the angular fields 9 ± and 9 2 
defined by 



Ax' 



1 



cos0 1 = -4- , sin^ = -Jc p if' 2 + c r </>' 2 , (65) 
& a' v 



and 



yj Cplfi' 2 + C r 1p' 2 ' 



sin 9 9 



Crip' 



yj Cplfi' 2 + C r 1p' 2 

(66) 

provided ip' and ip' do not vanish at same times. They 
define an instantaneous rotation matrix on the field- 
manifold 



M 



cos 9i sin 9i cos 9 2 sin 9 X sin 9 2 
- sin 9, cos 9 ± cos 9 2 cos 9 ± sin 9n 
cos( 



— sin 

— sin 



(67) 



for 9 1 and 9 2 , respectively. The rotated potential deriva- 
tives have been defined through 




= M 



(71) 



Note that since we are dealing with non-minimally cou- 
pled scalar fields, they are not the partial derivatives of 
the potential V with respect to the rotated fields [71| . 
however V a — L a V a remains the effective potential which 
sources the adiabatic field [see Eq. J57J)]. In the spher- 
ical basis, we recover explicitly that the entropy modes 
only source the comoving curvature perturbation since 
in Eq. (§21 one has ± c d V c 8T d = V l 8s l + V 2 Ss 2 . In order 
to compute the primordial power spectra and the cross- 
correlation for the different modes, it is convenient to 
recast the equations of motion for the original fields in 
terms of the rotated modes only [7(J [73, . 



which transforms the original fields into the adiabatic 
and entropy modes 




A5 X 
M I Cp/ 2 5tp 
cl^Sip 



(68) 



The evolution of the angular fields is readily obtained 
by differentiating Eqs. (|65|l and (|66|) and using the back- 
ground equations l|19f> to (|21|l 



a 2 V x 



+ a cos 9 1 I — — cos 



sin 9 2 



(69) 



sin Ma 



a 2 V 2 , 2 

; a COS I 

a' 



: sin 9 2 



— = COS Go , 

V* ) 

(70) 



3. Equations of motion 



The closed system of dynamical equations for the en- 
tropy and adiabatic modes can be obtained by express- 
ing the second order derivative of each mode with re- 
spect to the conformal time in terms of the others by 
means of Eqs. (O and (|ZU, using Eqs. lO to (J2U) as 
well as Eqs. (|69|) and (|7U[) and their derivatives. More- 
over, the adiabatic field 8a can be expressed in terms 
of the comoving curvature perturbations by means of 
Eq. (|60|l . which is a preferred observable for deriving 
the subsequent Cosmic Microwave Background (CMB) 
anisotropies (see Sect. HVf) . After straightforward but 
tremendous calculations one gets for the first entropy 
mode <5si 



10 



2a' ( a 2 V 2 a v . \ , 



x ( — — — sin 9 2 + cos 9 2 ) + sin 

7 /2 (Pw CQS 2 ^ + 



- A + a 2 Z tl 



11 + 4COS26*! + COS46*! - 
o ■ 2 „ a v ? 



- & 



- sin 2 9 2 



ha V a — — cos 2 H j=z 

Vc p c r / tan z X \ 



3 sin 2 6>! 
a 4 K 2 



Sin ^V ^ ^tan 2 ^ 



— — sin 9 2 + — -= cos 2 

. /Cp \/C r 



sin 9 2 



5s, + I 2a 2 V 12 + 



—9 cos 8 ± + cos 36>! 2 x 
2^ a V \ 

2a 2 V 2 f a 2 V t 3Ha' 



a 



sin( 



a 2 ^ 



— = cos 9 2 



tan 2 9 1 

'15- 8cos20 x +cos40! 2 

o a Vi 

4sin 2 6» 1 

2 - ' 



(72) 



2a 2 V 2 
sin X tan 8 1 \ 



a — cos 9 2 + — = sin 9 2 



ct' 2 Vtan 2 9 1 ^ tan6>! 



+ cr' 2 cos : 
where use has been made of 



a 12 cos 0, 
I i 

tan 2 6»! 
cos 29 



-| r. j sm20 2 H cos 20, 



Cr 



cr' 



- a 



I A 2q2V W 



'cr" \ 2a 2 
— - WJ <5a - — (K5 Sl + y 2 A S2 ) , 



a 2 F 2 



and, from Eq. I(63j) . 

a' 2 2W + a' 2 6Sl a- 
also introduced a rotated 2-form quantity (with respect to the field-manifold) from 



■5s 2 . 



We have 
to 



Vera 


v al 


v a2 


v 1<7 




v 12 


v 2a 


v 21 


v 22 



/ Ym. Ym v ^ \ 

A A^Cp A-^Cy 



A-yJ Cp Cp y j ) ■ , 

Vipx 



M 



-l 



(73) 
(74) 

the potential according 
(75) 



ensuring the symmetry properties 

\ 

Note that the special form of the effective potential in Eq. 1(11(1 has been 

_J_ L cos9 2 + sin8 2 ) (- sindM + cos^K) + (Ki - K.) 

cos 29 1 [ Vv^P v^r / 2 



= F 1<7 , V a2 = ^ 2ff , = 7 1B . 

used to perform some simplifications 



and 



V a2 = 



cos 6*i 



--^sin6» 2 + -^=cos6i 2 ) (- sin^Vi + cos^K) + sin 6^ 



In Eq. 172(1 . it was also convenient to introduce the exact potential derivative 

?%Wi6 2 ) {-sm8 l V l +cos9 1 V a ) . 



Zu = = Ki + cos 9 1 sin 9 t ( cos 9 2 + ■ 
OS! \ .fc^ 



(76) 



(77) 



(78) 



(79) 
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Similarly, for the other entropy mode one gets 

„ , 2er' / a 2 V 2 a,„ a,i, \ „ , 

5s'i + 2H5s' 2 - - — - — ^ + sin 9 2 - cos 9 2 5s[ 
tan6»! \ cr' 2 ^fc^ .frZ I 1 



^ 2 y 2a 2 V 2 ( a v 
22 tan 2 9 1 \ y/c^ 



— SI 



sm 2 



1 + sin 2 6i /2 { o-u, 
^ cr -sin( 



tan 2 0! 



a' 2 sin 2 6»! 

V 



(T COS 



?qA^£* S i n 2 # 2 + cos 2 # 2 - sin 2# ; 



— =: cos # 2 



o~s 2 



p^r 

'-3 + cos 20! 2 4 o 

-a K + 7" a ^ 

tan 1 



sin 2 6»! 



sm 6* 2 -E= cos 6> 2 + . - — — cos 6> 2 H sm 2 

y/Cr J sm 9 X tan t^j \ ^/c^ ^fc^ / 



(80) 



\ V f v / A 

2a 2 T/ 2 / a 2 ^ 3K(r'+a 2 K \ 
cr' 2 Uin 2 ^^ tan^ ) 



cr'' 



sin 0! tan X 



— - - sm 2# 2 H cos 26> 2 



2a 2 K w 



The equation of motion Eq. 1)63(1 for £ can also be expressed in terms of the entropy modes only by differentiation 



2H fa 2 V lx , a 2 V 2 A 2H 
— os, H -ds 2 



2W f tan 2^ , 2t ^ n . 
7 — (aVn-a 2 ^) 



sin^ [sin 2 l H ^— ) a 2 ^ + cos^ (cos 2 9 1 — - 

\ 008 201/ \ cos 26^ 

a 4 V 2 a 2 V, 



— — cos £/ 2 H 7= sm # 2 



COS (h 2t/ 

tan 1 
2H 



— sm 2 H -= cos tf 2 



•y/c^ / cr' 2 tan 0! cr 



/2 



a 2 ^ CT + cr' 



Us, 



(81) 



2 /x 

-< tan^a 2 !^ 



(-3 + COS20Q 2 2 

a K, 

2 8^26*! 



COS26*! - 7 2 / a v . a^, 

-a Kt I — sm9 2 ^cos^a 



a 2 F 2 / 2 



cr 



-/2 



a 2 ^ 



a 2 K, + f^+a 
tan 1 



nil s 



>Ss 2 . 



w,±± v 1 y ' u / J I 

In order to set quantum initial conditions during the inflationary eras, it is more convenient to use the Mukhanov 
fable [U Q — v'Q/7~L whose evolution is readily obtained from Eq. (|81|l 



vari 



Q" + 2HQ'+\ -A + a 2 V CTCT - 



— = sm 0., 



a 2 ^cos 2 ^ (— ^sin^ + ^Lcos^ 



cos 36>! - 5 cos 6»! 2 2 . 2 \ / a v 
a 2 V 1 + cos 2 Q x sm^i^Kj cos ^ 2 + 



— a 



(a 2 ^! - a 2 V aa ) 



sin ^ [sin 2 9 1 H ^— ) a 2 V t + cos 6» x (cos 2 9 X 

\ cos 2^! / \ cos 26*! 



a 2 V a 



cos 9 2 + —p= sin 9 2 



COS X 2 

tan( 



a 



a 2 V 2 ^sin6» 2 + — ^=cosf 



- 2 1 tan 9 1 a 2 V 12 



cr' 2 tan 9 1 cr' 2 
-a 2 V a 



a 2 V a + cr' Un-^ 



(-i + cos26 1 ) 2 2 cos20!-7 2 1 / a v . ou, 

— „ . n „ — —a^ H a 2 V a sm9 2 ^ cos^ 

2 \ ,/&T ^/c7 



2 sin 26>! 



(82) 



a 2 y 2 / , a 2 V, 

' < ~ 2 V„ + 1 : - + a '[m- — 



t'2 



,2t/ / W r 

tan&i 



- cr 



Us 2 - 



The equations J75J, (|80|l and l|82|l form a closed system moving curvature perturbation ( in Eqs. I|72|l and 1801) 
of equations. Indeed the source terms involving the co- 
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can be expressed in terms of Q as 



H s 



2 V 



Q' 



a 2 V a 



H 



(83) 

One can check that Eqs. lf721) . (|8*Tj) and match with 
the ones obtained for two fields in Refs. |46l 1771 ] for Sip — 
■0 = 0, whereas Eq. (|80[) decouples and admits Ss 2 = as 
a solution, up to transformation 9 — > 7r/2— 9 t . As shown 
in Sect. [HJ since the dilaton ip relaxes toward i/i = 
rapidly, we expect to recover a two field like behavior 
also at the perturbation level. In the next sections we 
numerically solve the system of equations l)72|) , l|8U|l and 
(|S2|) to compute the resulting power spectra at the end 
of inflation. 



4- Numerical results 

For a given background evolution (see Sect. [HI, and 
in the Fourier space with respect to the comoving spa- 
tial coordinates, the solutions of Eqs. (|72|) . (|%T|) and 
are uniquely determined by the initial values of 5s u Ss 2 , 
Q and their first order time derivative. Motivated by 
the quantum origin picture of the cosmological perturba- 
tions, in the Einstein Frame the scalar and tensor modes 
are decoupled and we can treat 8s u Ss 2 and Q as inde- 
pendent stochastic variables deep inside the Hubble ra- 
dius HHH Sill. Following the method of Ref. 
since we are interested in the two-point correlation func- 
tion between the different rotated fields <5s 1; Ss 2 and Q (or 
£), we will consider only the solutions obtained by start- 
ing with one mode in the Bunch-Davies vacuum while the 
others vanish, and the other mode permutations. In addi- 
tion to providing enough information to compute all the 
power spectra, they clearly exhibit the sourcing effects 
of one field from the others. For numerical convenience, 
each Fourier mode k will be assumed to appear at a given 
time rj q such that 



a 



(84) 



where k is the comoving wave number and C q a con- 
stant verify ing C n 3> 1 and characterizing the decou- 
pling limit |46l l96j . For a Bunch-Davies vacuum, the 
value of C q does not usually change the power spectra, 
provided it is big enough to allow free wave solution 1 . 
For [i s = a^/2k z / 2 Q (and respectively fi s — aV^k^/^Ssi, 
Us = aV2k 3 / 2 5s 2 ), we therefore set at 77 = r] q |l00l llOlj 



k 



(85) 



1 Th is is no longer true for trans-Planckian initial conditions 



In Fig. the numerical solutions of Eqs. (|72|) to l|82|) 

have been plotted for the same background evolution as 
in Fig. [21 and for three different values of the comoving 
wave number k with k x < k 2 < k 3 . In the top frame, the 
evolution of the Hubble parameter and the three phys- 
ical wavenumbers k/a are plotted as a function of the 
number of efold before the end of inflation ("bfold" in 
the following). A subtlety occurs however for the "very 
large" wavelength modes due to the existence of a non- 
inflationary era when the background evolution is domi- 
nated by the rolling of the ip field (see Sect. Inland Fig.[2J). 
Indeed, if the initial conditions in the five-dimensional 
setup are such that the initial value of ip is far from van- 
ishing, then there are some perturbation modes which 
initially are super-Hubble and for which setting quantum 
initial conditions does not make sense. However, since 
we are interested in enough efolding to solve the flat- 
ness and homogeneity problems, such large wavelength 
modes are generically still super-Hubble today and thus 
non-observable. On the other hand, one cannot exclude 
that for a small number of efold, let's say 60 for instance, 
some of the large wavelength perturbation modes enter- 
ing the Hubble radius today may have been created dur- 
ing the transition period between the ip dominated phase 
and the inflationary one, as k 1 in Fig. [3] Such a mode 
starts initially under the Hubble radius, but there is a 
maximum value of C q for which Eq. I|84|l has not solu- 
tion (because ry q would be too small) . Our prescription is 
to discard the models for which observable perturbation 
modes today would have been created on Hubble, super- 
Hubble and sub-Hubble scales with C q < 100 initially. 
From the second to the bottom frame in Fig. |3J the evo- 
lution of |(5(7i-?i ond )|, \Ss 1 (n~n cnd )\ and |<5s 2 (n-7i en d)| 
are plotted, respectively. The straight, dashed and dot 
lines correspond to the three independent initial condi- 
tions for the modes, Q in a Bunch-Davies vacuum and 
5s 1 = 5s[ = Ss 2 = 5s' 2 = 0, Ss 1 in a Bunch-Davies vac- 
uum and Q — Q' — 5s 2 = Ss 2 = 0, . . . , respectively. 
Moreover, in each frame we have plotted the three wave- 
length modes k 11 k 2 and k 3 from the left to the right. 
In Fig. El we have plotted the real and imaginary parts 
of the comoving curvature perturbation ( = QH/cr' for 
the mode k 1 associated with the Q solution of Fig. [S] 
and the background evolution of Fig. [21 This plot shows 
the transition between the oscillatory behavior when the 
mode is under the Hubble radius and its damping on 
super-Hubble scales. However, since we are precisely in 
a multifields case, even on super-Hubble scales the co- 
moving curvature perturbation slightly evolves sourced 
by the entropy modes till the end of inflation. Also, the 
rapid change in behavior of the Hubble parameter during 
the transition between the ip dominated era and the sub- 
sequent inflationary eras induces also a change of slope 
in the evolution of all the perturbation modes (see Fig. O 
and Fig. El around n — n on( j — —70). This can be under- 
stood by first noting that before the transition, as long as 
e > 1, the physical wavelength of the mode goes deeper 
into the Hubble radius with time whereas as soon as e < 1 
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FIG. 6: Real, imaginary and modulus of the comoving cur- 
vature perturbation mode ( = QTL/a' — Q/a corresponding 
to the fei wavenumber, and model parameters, of Fig. El The 
mode behaves as a free wave under the Hubble radius whereas 
it is damped after Hubble exit while sourced by the entropy 
perturbations. Note also the change of slope of the modulus 
during the background transition between the ip dominated 
expansion and the subsequent inflationary era (see Sect.lITl. 




FIG. 5: Evolution of the adiabatic and entropy perturbations 
in bfold time (efold before the end of inflation), for three 
wavenumbers k ± < k 2 < k 3 . In the top frame, the Hubble pa- 
rameter and the physical wavenumbers k/a are plotted. The 
next three frames shows the evolution of each Q, 5si and 
8s 2 modes, respectively, obtained from the three independent 
initial conditions where only one of the mode is in a Bunch- 
Davies vacuum and the others vanish. Only the wavenumber 
leaving the Hubble radius close to the background 4> domi- 
nated expansion exhibits strong couplings between the adia- 
batic and entropy modes. The background model parameters 
are the same as in Fig. Inland we have chosen fh c = 10 -6 . 



(see Fig. EJ). Note that in Eqs. to (JHOJ, the Hub- 
ble parameter appears also in the effective potential of 
the modes. Its rapid (exponential) evolution during the 
relaxation of ip enhances the coupling between the en- 
tropy and adiabatic modes. Also, at the transition, the 
background fields are no longer in a stationary regime as 
discussed in Sect.|Hl and especially one expects 9[ and 9' 2 
to increase and reinforces the coupling between the adi- 
abatic and entropy modes. This effect is similar to the 
one observed for double-inflation models 0]. For the 
modes which evolve out of this stage, the ip field and its 
perturbations are negligible, i.e. 9 2 ~ and the second 
entropy mode Ss 2 decouples and decreases exponentially 
[see Eqs. JHHJ, (JHOJ an d the bottom frame in Fig.0. The 
dynamics is close to a two fields like regime where the 
friction dominated evolution leads to a weakly coupled 
entropy mode Ss 1 0, |97|. As can been seen in Fig. [51 
after Hubble exit, for the wavenumbers k 2 and fc 3 , the 
entropy perturbations Ss 1 are quasi-frozen as it would be 
for a light scalar field. 



the wavelength becomes closer and leaves the Hubble ra- 
dius. Moreover the amplitude of the sourcing effect of 
the entropy modes on the comoving curvature perturba- 
tion is enhanced. As can be seen in Fig. and Fig. El the 
adiabatic mode is mainly produced by the entropy modes 
for k = k ± whereas this effect progressively disappears for 
the bigger wavenumbers k 2 and k 3 . This comes from the 
fact that only the k± mode evolves and leaves the Hubble 
radius in a regime where the ip field is not yet vanishing 



C. Tensor modes 

In the Einstein frame, the scalar and tensor degrees of 
freedom are decoupled. Therefore the equation of evolu- 
tion for the tensor modes remains the same as in General 
Relativity. For a flat perturbed FLRW metric 

d.s 2 = -a 2 dry 2 + a 2 (<% + h tj ) dx'dx 3 , (86) 
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FIG. 7: Real, imaginary and modulus of the tensor perturba- 
tion mode corresponding to the k 2 wavenumber. The Hubble 
parameter is also shown and we recover the standard result 
\fl T /a\ — Hk with Ht its value at Hubble exit. Note however 
that it takes few more efolds for the mode to remain constant 
which means that this approximation is valid only for a slow 
varying Hubble parameter (as during the slow-roll stage). 



where hij is a traceless and divergenceless tensor 

S^hij = 5 tk d k h t] = 0, (87) 

one gets [H [To2 

h'lj + 2Hh' tJ - Ahij = 0. (88) 

This equation can be numerically solved for each po- 
larization state h by setting the initial quantum tensor 
modes /i T = k 3 ^ 2 ah in the Bunch-Davies vacuum [see 
Eq. H85fl ]. hi Fig. [7| their evolution has been plotted 
for wavenumber k 2 together with the evolution of the 
Hubble parameter. As expected for a massless field, the 
gravitational waves freeze on super-Hubble scales with 
|/x T /a| ~ Hk where Hu is the Hubble parameter at Hub- 
ble exit. Note however that the freezing occurs a few 
efolds after Hubble exit. 



D. Primordial scalar and tensor power spectra 

From the scalar and tensor modes evolution obtained 
in the previous sections, the primordial power spectra are 
readily obtained from the values taken by the adiabatic 
and entropy perturbations at the end of inflation. In 
the Fourier space and Einstein frame, the scalar power 
spectra are computed according to the method used in 
Ref. Htf|, i.e. 

^ b = |^EK(fc)]*[^(fc)], (89) 

rn 

where the observable perturbation modes v a stands for 
C, Ss 1 /a and Ss 2 / '<r, whereas the index "m" refers to the 



three independent quantum initial conditions. The ten- 
sor power spectrum reads 

o L.3 

Vh (k) = ^ T \h(k)\\ (90) 

where the polarization degrees of freedom have been in- 
cluded. 

In Fig. |H1 we have plotted the typical power spectra 
associated with the solutions computed in Sect. IIII B 41 
The adiabatic power spectra is dominant compared 
to the two entropy ones Vs t and Vs 2 - Furthermore, the 
second entropy power spectrum Vs 2 is strongly damped 
and blue tilted compared to Vs 1 ■ In fact, this behavior is 
expected since the second entropy mode Ss 2 follows the 
tjj evolution during inflation and vanishes exponentially 
with the total number of efolds. As a result, the earlier 
the Ss 2 modes cross the Hubble radius, the longer they 
sustain an exponential decay. Their power spectrum at 
the end of inflation is blue tilted since the smaller scales 
are accordingly less damped. However, since the observ- 
able perturbation modes are expected to cross the Hub- 
ble radius around 60 efolds before the end of inflation, 
the entropy power spectrum Vs 2 cannot be significant at 
the end of inflation. In particular, the expected cross- 
correlations for fcj-like modes discussed in Sect. IIII B 41 
are also washed out by such a damping. They can still 
be seen in Fig.|Hlas a bump in the cross-correlation power 
spectra V^s 2 and Vs 1 s 2 : but with an extremely small am- 
plitude. It is important to recall that the domination of 
tp during inflation has to occur before the last 60 efolds 
of inflation since it drives a non-accelerated expansion 
in our model. As seen in Sect. |Hj this comes from the 
fact that c r + c p = 6 whereas in a more generic scalar- 
tensor model with c r + c p < 9/2 the transition from 
-0-dominated to (/3-dominated inflation could be in the 
observable range and may generate a si gnifi c ant powe r 
spectrum for the second entropy mode |l03t Il04 Il05| . 
Nevertheless, we will accordingly be focused in the fol- 
lowing on the present boundary inflation model where 
only the adiabatic, tensor and first entropy modes are of 
interest for the CMB j4(|. 

The hierarchy in the power spectra amplitudes can 
be qualitatively understood by noting that the entropy 
mode i5si behaves almost like a free massless field for 
moderate values of the coupling constant c p . Indeed, as 
can be seen in Fig. 5s ± is weakly coupled to the adi- 
abatic mode and remains almost frozen once it becomes 
super-Hubble. Neglecting the sourcing effects from the 
adiabatic modes, as well as the variation of the Hubble 
parameter during the few efolds after Hubble exit, one 
gets at the end of inflation 

fc 3/2 |^le„ d -^ (91) 

where Hf. is the dimensionless Hubble parameter kH 
at Hubble exit for each mode of wavenumber k. From 
Eqs. i|39[l and Q89p. keeping in mind that e ~ 1 at the 
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FIG. 8: Typical power spectra at the end of inflation for the 
perturbation modes of Fig.|^|and Fig.[7|which becomes super- 
Hubble after efolds. The entropy power spectra are sub- 
dominant compared to the adiabatic one. 



end of inflation, the resulting power spectrum varies as 



1 

16' 



(92) 



Similarly, assuming that the adiabatic mode Q gets mass- 
less field-like fluctuations at Hubble exit, in the weakly 
coupled regime the comoving curvature perturbation re- 
mains constant afterwards and the adiabatic power spec- 
trum can be approximated by 



8^ 2 e k ' 



(93) 



where e k is the first slow- roll parameter at Hubble exit. 

One can check in Fig. |H| that Eqs. J32J) and pro- 
vide a good approximation of the relative power spectra 
amplitudes. Moreover, these equations show that the 
generated power spectra do not depend significantly on 
the initial values of the background fields which drive in- 
flation. Indeed, as thoroughly discussed in Sect. [H] the 
fields are rapidly attracted toward a friction dominated 
regime during which their evolution does no longer de- 
pend on their initial value. Obviously, the initial condi- 
tions determine the total number of efolds, but at a given 
cfold before the end of inflation one may not expect to 
be strongly influenced [see Eq. (53} ] ■ However, the ampli- 
tude of the perturbations at Hubble exit directly depends 
on the Hubble parameter value which involves both the 
bare mass m c of the matter field and the conformal fac- 
tor A. In order to disentangle their respective effects, it 
is more convenient to consider the effective mass of the 
matter field at the end of inflation 

TOeff = A cnd m c = A cnd nm c . (94) 

In Fig. |31 the slow-roll and Hubble parameters have been 
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FIG. 9: Evolution of the first slow-roll parameter e and the 
dimensionless Hubble parameter H for different initial values 
of the background fields. The <p coupling constant has been 



set to c D 



0.1 and the effective mass of the matter field at 



the end of inflation is m e ff = 6.2 x 10 . The existence of a 
friction dominated regime for the background fields leads to 
observable power spectra which are insensitive to the initial 
values of the background fields. 



plotted for several arbitrarily chosen value of Aj n iXini and 
i/'ini, the other parameters c p and m e fj being fixed. As the 
plots show, the last efolds of evolution cannot be differen- 
tiated and we have checked by a direct computation that 
this is also the case for the corresponding power spectra 
at the end of inflation. In fact, the only effects ^4iniXini 
and ^inj might induce concern the transition at the end 
of the i/j-dominated expansion, which is, as explained be- 
fore, hardly observable for the CMB. 

The remaining degrees of freedom in the power spec- 
tra are the effective mass rh e g and the moduli coupling 
constant c p . These parameters are expected to have sig- 
nificant observable effects. On one hand, the effective 
mass fixes the overall amplitude of the Hubble parame- 
ter and is therefore directly related to the amplitude of 
the primordial perturbations. On the other hand, the 
coupling constant c p encodes how much the conformal 
factor may run during the generation of the observable 
perturbations. From Eq. i|ll[l . such a running renders 
the potential steeper and one may expect bigger tilts for 
the power spectra. 

In Fig. ^3 the slow-roll and Hubble parameters have 
been plotted for several values of c p , the effective mass at 
the end of inflation being fixed. The slopes of these two 
functions are effectively steeper for larger values of the 
coupling constant and the resulting power spectra may 
be qualitatively guessed from l|92|l and (|93|1 . The direct 
computations of the power spectra at the end of inflation 
have been plotted in Fig. II II and confirm this behavior. 

In order to dress a qualitative understanding of the 
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FIG. 10: Evolution of the first slow-roll parameter e and the 
dimensionless Hubble parameter H for c p = 0.1, c p = 0.01 



and 



= 0.001. The effective matter field mass is m e s — 



6.2 x 10 . The running of the conformal factor during the 
generation of the observable perturbations leads to more tilted 
power spectra (see Fig. 111! . 
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FIG. 12: The adiabatic power spectrum at the end of a weakly 
coupled boundary inflation compared to the one obtained in 
single field chaotic inflation. The dotted line is the first order 
slow-roll approximation of the chaotic model spectrum around 
the pivot scale fc* . The deviations between the weakly coupled 
boundary and chaotic models come from the slight running 
of the conformal factor. 



trum simplifies to [10 
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FIG. 11: Dependence of the power spectra at the end of infla- 
tion with respect to the moduli coupling constant c p . Larger 
slopes are obtained for higher values of the coupling constant, 
as expected from the running of the conformal factor during 
the generation of the perturbation modes (see also Fig. I1UI . 



above-mentioned effects, we have compared in Fig. IT21 
the adiabatic power spectrum generated during a weakly 
coupled boundary inflation with c p = 3 x 10~ 3 with the 
one associated with the standard single field chaotic infla- 
tion (c p = 0). Moreover, we have plotted the first order 
slow-roll approximation of the chaotic model around a 
pivot scale fc*. In this limit, the adiabatic power spec- 



where all the parameters are evaluated at the pivot scale. 
The constant C sr ~ —0.73 and e 2 = dlne/dn is the 
second slow-roll parameter. The deviations observed in 
Fig. El between the slow-roll and single field power spec- 
tra come from the natural running of the spectral index 
which is not grabbed by the first order slow-roll approxi- 
mation. The differences between the chaotic and bound- 
ary spectra are due to the slight running of the conformal 
factor obtained for c p = 3x 10 -3 . More qualitatively, for 
non- vanishing values of the moduli coupling constant, the 
effective mass of the matter field is not constant during 
the generation of the observable perturbations and in- 
duces deviations with respect to a single field model. 

As previously pointed, Eqs. I|92|l and l|93|) are relevant 
only if the modes do not evolve significantly after Hubble 
exit. This is obviously the case for the single field model, 
but certainly no longer true for high values of the moduli 
coupling constant. In Fig.^J we have plotted the ratio 
between the computed power spectra for several values 
of c p . Some deviations show up for c p > 0. 01 in the r atio 
tensor to scalar and in the entropy modes |1071 Il08| . 

In the next section, after setting up a toy cosmological 
model, we use the above numerical method to derive the 
resulting CMB anisotropics and perform a Monte Carlo 
Markov Chains exploration of the model parameter space 
given the first year WMAP data. As expected from the 
behaviour of the primordial power spectra, we find both 
the effective mass of the matter field at the end of in- 
flation and the moduli coupling constant c p to be con- 
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FIG. 13: Ratios between the power spectra of Fig. 1111 and 
their evolution with respect to the moduli coupling constant 
c p . The "consistency check" Vh/Vc, = 16e is recovered in the 
weak coupling limit while significant deviations appear for 
c p = 0.1 due to the sourcing effects between the modes after 
Hubble exit. 



strained. 



IV. CMB ANISOTROPIES 

A. Cosmological framework 

As mentioned in Sect, [n] in the brane world picture, 
the moduli ip and if are associated with the position of 
the branes in the five-dimensional spacetime while the 
field x lives on the brane supposed to be our universe. 
As a result, a natural setup is to consider that the moduli 
fields remain in the late-time cosmology and in particu- 
lar today whereas the radiation, matter and dark matter 
sectors are sourced by the decay of the field \ through a 
reheating period [H El EH (see Fig. EJ) . In the follow- 
ing we consider a toy cosmological background model by 
assuming instantaneous reheating [109| . In this respect, 
the energy density at the beginning of the radiation era 
is completely sourced by the energy density of the matter 
field x a t the end of the inflationary period. This allows 
an estimation of the scale factor at the end of inflation 

m ^~_ m ^ + I m -^L ; (9 6 ) 
a a cq 4 p cn d 

where the index zero refers to today, "eq" at the equiva- 
lence between radiation and matter, and with the energy 
density p ond = p x . 

On the other hand, since we are dealing with a 
multiscalar-tensor theory in the late-time cosmology, and 
especially today, the evolution of the conformal factor 
and its first and second gradients are constrained by var- 
ious solar-system, astrophysical and cosm olog ical mea- 
surements H3, Ea E3, EE EE E3, E3, E3, E3| a s 



discussed in Sect. [H] the field if) is rapidly driven toward 
zero and we will safely assume in the following that it 
is indeed the case during the radiation and matter eras. 
From Eq. J7J, the only non- vanishing first conformal gra- 
dient is a v , which remains constant during the cosmolog- 
ical evolution. The strongest constraint comes from the 
solar system m l62 |- Fr om the Cassini spacecraft mea- 
surements, one gets |llOf 

e ab a a a b < 5 x 1(T 7 => c p < 2 x 1CT 5 . (97) 

solar solar 

Note that once i/j = 0, the only non- vanishing second con- 
formal gradient is poorly constrained. This comes 
from the fact that the post-Newtonian parameter 7 pn in- 
volves only the combination a a a b /3nh which vanishes in 
that case, whatever the value of c r [5(J]. The other con- 
straints coming from the v ariation of the "Cavendish" 
gravitational constant |lll| 

/&, = K 2 A 2 {l + 2£ ab a a a b ), (98) 

are found to be satisfied once Eq. ij!T7|) is. The cosmo- 
logical constraints on scalar-tensor gravity are less strin- 
gent than the solar-system ones but allow one to probe 
the cosmic times. T hey lead to the two-sigma upper 

bound m m in ei ei m m m m ei 

e ab a a a b < 2 x 1CT 3 => c p < 8 x 1(T 2 . (99) 

cosmo cosmo 

In the following, we are interested in the observable 
consequences the previously discussed boundary infla- 
tion eras and their resulting primordial power spectra 
may have on the CMB. The fact that the x ne ld de- 
cays into radiation, matter and dark matter implies that 
there are no observable entropy modes between the pro- 
duced cosmological fluids. Of course they are present 
between the cosmological fluids and the moduli fields 
ip and ip, however we will assume that any back reac- 
tion effects on the minimally coupled fluids can be ne- 
glected as soon as the theory does not deviate too much 
from General Relativity after inflation. As a result, only 
the adiabatic and tensor primordial power spectra de- 
rived in Sect. ITTTI source the observed cosmological per- 
turbations. The current cosmological data being sensi- 
tive to the tilt and amplitude of th e primordial power 
spectra [HI Ell El Ell Ell EH , one may expect to 
probe the model parameters involved. In particular, since 
for c p = our model matches with single field chaotic in- 
flation, this parameter will be used to analyse whether 
the boundary inflation model is favored or excluded by 
the data compared to a single field model. 

Keeping in mind that the present model has to verify 
Eq. (|97|l . we will nevertheless only use a weak late-time 
cosmology upper bound Cn < 0.2 as a prior in the follow- 
ing CMB computations El E3 • Indeed, previous deriva- 
tions of the CMB anisotropies in scalar-tensor gravity 
theory have been focused on the late-time modifications 
the non-minimally coupled scalar fields may produce by 
assuming standard power law primordial power spectra. 
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FIG. 14: Correspondence between the wavenumber today k 
(Mpc -1 ) and the bfold of Hubble exit. The background model 
is the same as in Fig. |H1 




FIG. 15: Dependence of the TT power spectra with respect 
to the moduli coupling constant c p , at fixed effective matter 
field mass m e ff = AmdW- c . 



The present analysis being precisely concerned with the 
primordial stages and the generation of the cosmological 
perturbations, it may give a complementary view of the 
effect expected on the CMB in scalar-tensor theories. 

Under these assumptions, one can approximate K e g = 
kA ~ 

^cav and Eq. I|96[) can be further simplified 
into [n|| 

ln^c,hn(J^K eS H )-hn^^, (100) 
a B 2 \ v J 4 3 

where a , r2 ra( j and H are, respectively, the scale factor, 
the total density param eter of radiation and the Hubble 
parameter today |l02j |. Note that n^ s p cn d ~ n 4 p C nd in 
the Einstein frame for small variations of the conformal 
factor during the late-time cosmology. From Eq. (|100fl . 
we can relate the wavenumber of the observed perturba- 
tions today to the corresponding comoving wavenumber 
during their primordial generation 

£ = C c |^e-(— >, (101) 
rt ti a en d 

where k is the comoving wavenumber today in units 
of Mpc~\ H — kH and C c a conversion unit constant 
C c ~ exp(— 130). In Fig. ^] we have plotted the corre- 
spondence between k and the bfold of Hubble exit during 
inflation for the same model parameters as in Fig. El an d 
with a = 1. 

B. CMB power spectra 

In order to compute the multipole moments of the 
CMB anisotropies we have used a modified version of 
the CAMB code [U. From Sect. |H| and Sect. EH 
in addition to the usual cosmological parameters, the 



CMB anisotropies in the boundary inflation model un- 
der scrutiny are characterized by A- m i, ip- m i, and c p for 
the non-minimally coupled sector as well as Xini and fh c 
in the matter sector [see Eqs. © and ©]■ However the 
absolute value of the conformal factor A is non-observable 
and can be arbitrarily set at a given time, e.g. Aj n i = 1. 
As a result, it is more convenient to parametrize the 
boundary inflation eras in terms of the rescaled param- 
eters: ^liniXini, m Q g ee A cnd m c = K cS m c , with ip ini and 
c p unchanged. For given values of these four parame- 
ters the background solution is computed and used to 
setup the cosmological framework by means of Eqs. I|100|) 
and (|101f> . Then, each observable wavenumber required 
in the primordial adiabatic and tensor power spectra is 
computed from its corresponding quantum initial state 
as presented in Sect. Hill and thereby used to derive the 
resulting CMB power spectra. According to the discus- 
sion in Sect. IIIIDl we do not expect observable effect on 
the CMB coming from ^4i n iXmi and ipini since for a wide 
range of values they do not modify the primordial power 
spectra. However, these parameters have been kept free 
in the following in order not to reject the rare cases for 
which the end of the -^-dominated expansion era would 
precisely occur during the generation of the largest ob- 
servable perturbation modes. On the other hand, Fig. 1151 
confirms that varying the moduli coupling constant c p is 
not innocuous for the temperature angular power spec- 
trum. However, as previously mentioned, the effective 
matter field mass rh c g is also directly involved in the am- 
plitude of the primordial power spectra through Eq. I|2rt(l 
and one may expect a degeneracy between these parame- 
ters on their respective CMB influence. In addition, both 
c p and m e ff have been shown to modify the tilts of the 
primordial power spectra in Sect. iHIDl 

In the next section, a Monte Carlo Markov Chain 
(MCMC) exploration of the cosmological and primor- 
dial parameter space if performed by using the first year 



19 



WMAP data. Note that the method we are using relies 
on a full numerical scheme which does not involve any 
approximation of the primordial power spectra. 

C. MCMC exploration 

Following the method of Ref. [70|. we consider a pa- 
rameter space involving a minimal set of cosmological 
parameters: the ratio of the sound horizon to the angu- 
lar diameter distance 9 (related to the reduced Hubble 
parameter h), the density parameter of baryons f2b and 
cold dark matter f2 c , the optical depth r (or the redshift 
of reionization z rc ); as well as our primordial parameters: 
AniXini, ipini> c p an d w e ff. Note that the four primordial 
parameters fix the overall and relative amplitudes of the 
primordial scalar and tensor power spectra. The MCMC 
computations have been done by using the COSMOMC 
code |70j calling the modified CAMB version based on our 
inflationary code, and given the first year WMAP data 
and the associated likelihood code [6{| [6(1 [62l . 

In order to check the relevance of our full numerical ap- 
proach, we have first performed a MCMC exploration on 
the single field chaotic inflation model which is obtained 
by fixing -0; n i = c p = 0. Indeed, the current constraints 
on the cosmological parameters using the WMAP data 
usually assume either power law or slow-roll a pproxi- 
mated primordial power spectra |l!3t Ill4l Ill5| . Since 
our method goes further it may be interesting to check 
its consistency with the current existing bounds. More- 
over, the single field chaotic model will be our reference 
model to discuss the more complex features associated 
with the boundary inflation model. 



1. Chaotic model 

We have used standard prior distributions for the base 
cosmological parameters (see Ref. |70|) whereas wide top 
hat priors have been chosen for the chaotic model pa- 
rameters: 20 < Xini < 1000, and -10 < log(m c ) < 0. 
The lower limit on the initial matter field value is set in 
order to get the right order or magnitude of the mini- 
mal total number of efolds required to solve the flatness 
and homogeneity problem (~ 60). In addition, in the 
COSMOMC code, we have coded a "hard prior" rejecting 
any model for which observable perturbations cannot be 
initially set in a Bunch-Davies vacuum with C q = 100: 
for the chaotic model this can occur when there is not 
enough efolds of inflation [see Eq. 1)84(1] . The higher limit 
has been chosen in order to avoid prohibitive computa- 
tion time which occurs when the total number of efolding 
become very large. The prior on log(m c ) is chosen to con- 
tain the value required to get the right amplitude of the 
CMB anisotropics. The obtained posterior probability 
distributions for the base and derived cosmological, as 
well as the primordial parameters are plotted in Fig. 1161 
and Fig. El They correspond to 50000 samples for which 
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FIG. 16: The ID marginalized posterior probability distribu- 
tions in the pure chaotic inflation model. The dotted lines are 
the associated ID mean likelihoods. 
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FIG. 17: la and 2a isocontours of some correlated 2D 
marginalized posterior distributions in the pure chaotic in- 
flation model. The colormap traces the 2D mean likelihoods. 
The inflaton mass fixes the amplitude of the primordial per- 
turbations and is constrained around Grand Unification val- 
ues. 



the errors on their shape do not exceed 3%. Note that 
they do not rely on any slow-roll approximation but only 
on the linear perturbation theory and the cosmological 
setup of Sect. HvH 

Firstly, the constraints on the base cosmological pa- 
rameters are found to be cons istent with the current state 
of the art [111 El [ill EHI • For the chaotic model pa- 
rameters, the inflaton mass m c is well constrained, as 
expected from a parameter involved in the normalization 
of the primordial power spectra, and we find, at 2a level 
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-5.24 < log(Km c ) < -5.18. (102) 

As expected, the initial value of the field is found to 
be unconstrained. Concerning the overall ability of the 
chaotic model to fit the data, the best fit is obtained 
with a likelihood of — ln(£) ~ 714.4 for 1343 degrees of 
freedom which render the chaotic model slightly more 
favored than its slow-roll approximated version. 

In the next section, a MCMC exploration is performed 
for the boundary inflation model. Since in the limit c p — > 
this model matches with the chaotic one, the posterior 
probability distribution of the coupling constant c p traces 
how favored or disfavored the boundary inflation model 
is with respect to the single field chaotic model. 

2. Boundary inflation model 

The full set of primordial parameters is now consid- 
ered: Ai n iXini, V'ini, c p and m e ff- Compared to the chaotic 
model, there are two additional parameters ipi n i and c p , 
and a rescaled one: m c g = A cn dKm c . Remember that the 
latter no longer refers to the bare mass of the matter field 
but to its effective mass at the end of inflation. A flat 
prior has been chosen for the moduli coupling constant: 

— 10 < log(cp) < —0.75. The lower limit corresponds to 
the decoupling between the ip and x fields for which the 
boundary model is equivalent to single field chaotic in- 
flation, whereas the upper limit correspond to the weak 
post-inflationary bound above which the moduli influence 
the late-time cosmology (see Sect. IIV . Concerning 
the ipi n i values, we have set a conservative top hat prior 
< ipini < 1 since, according to our C q = 100 prescrip- 
tion, larger values of ip cannot be observed from the CMB 
point of view (see Sect. HUB 40 . The other priors for the 
base cosmological parameters and the effective mass are 
the same as for the chaotic model. Their correspond- 
ing posterior distributions have been obtained for 160000 
samples and are plotted in Fig. 1181 and Fig. 1191 where the 
marginalized distributions obtained for the chaotic model 
are also superimposed (red dashed lines). 

On one hand, the robustness of the cosmological pa- 
rameters with respect to the primordial power spectra is 
recovered: there are no significant deviations between 
their posterior probability distributions in the chaotic 
model and in the boundary inflation model. 

On the other hand, we did not find a better fit of 
the WMAP data with the boundary inflation primordial 
power spectra. More precisely, the best fit has the same 
likelihood than the one obtained with the chaotic model: 

— ln(£) = 714.4 but now for 1341 degrees of freedom. 
Statistically, the boundary inflation model is disfavored 
by its additional two param eters which do not help to im- 
prove the fit [7(J, ll2Ct Il2l| . As expected from the previ- 
ous discussions, the initial values of the background fields 
AniXini and ^ini have not observable effect and are con- 
sequently unconstrained by the data. Furthermore, the 
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FIG. 18: The ID marginalized posterior probability distribu- 
tions (black solid lines) and the ID mean likelihoods (dotted 
lines) in the boundary inflation model. Only the base cosmo- 
logical and primordial parameters are represented, together 
with the chaotic model posteriors of Fig. llfiK 'red dashed line). 
These posteriors have been derived under the prior choice 
log(cp) < 0.2 for which the moduli do not strongly influence 
the late-time cosmology. 
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FIG. 19: la and 2a isocontours of the correlated 2D marginal- 
ized posterior distributions in the boundary inflation model. 
High values of the coupling constant c p lead to bigger tilt in 
the primordial power spectra and are strongly cut. The de- 
generacy between the effective mass of the matter field and 
the moduli coupling constant is seen on the right plot. 



marginalised probability distribution of the coupling con- 
stant c p in Fig.ll8lexhibits a plateau for c p < 10~ 3 . Since 
in the decoupling limit c p — * the boundary model iden- 
tifies with the single field chaotic model, the existence 
of such a plateau shows that the WMAP data are not 
sensitive to the deviation induced by the moduli during 
inflation as long as c p < 10~ 3 . As a result, from the data 
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point of view, all these inflationary models, boundary 
and chaotic, cannot being distinguished in that regime. 

However, the rapid decay of the c p marginalised prob- 
ability distribution for values larger than 10~ 2 signals 
the sensitivity of the data in this regime, which disfavour 
the corresponding boundary inflation models. More pre- 
cisely, one obtains the two-sigma level upper bound for 
the moduli coupling constant 

log(cp) < -1.64, (103) 

which can be recast to an upper limit for the scalar-tensor 
conformal gradient 

a 2 = e ab a a a b < 6 x 1(T 4 , (104) 

at 95% confidence level. In terms of the post-Newtonian 
parameters more commonly used to constrain scalar- 
tensor theories jH^], we have 

4a 2 

1_7= TT2a^ <2Xl ° (1 ° 5) 

In fact, this constraint comes from the effects discussed 
in Sect. IIII Dl the tilt of the primordial power spectra 
in the boundary model increases with the value of the 
the moduli coupling constant due to the running of the 
effective matter field mass. Since the WMAP data are 
compatible with flat primordial power spectra, it is there- 
fore not surprising that they end up constraining devia- 
tions from scale-invariance. This effect also explains why 
the marginalised probability distribution of the effective 
mass m e ff = 4 cn( im c is skewed to lower values. For a 
given value of m e g, larger values of c p implies a more ef- 
ficient running of the conformal factor A during inflation, 
and in particular larger values of Afh c at the time when 
the observable perturbations have left the Hubble radius 
(see Fig. ITTI and Fig. IT2"|l . This degeneracy between c p 
and m e g is clearly seen in Fig. 1191 the two-dimensional 
marginal probability is distorted in such a way that the 
overall amplitude of the primordial perturbations fit with 
the one measured by WMAP. 

Let us emphasize again that Eq. \ 104(1 holds at the 
time of inflation and comes only from the shape of the 
primordial power spectra. In fact, for a given matter 
sector, this bound certainly applies to any inflationary 
scalar-tensor models involving a running conformal fac- 
tor during the generation of primordial perturbations. 
On the contrary, scalar-tensor models exhibiting a relax- 
ation mechanism toward General Relativity are expected 
to evade this bound: as it is the case for the moduli ip 
in our setup, such scalar gravity fields would reach the 
minimum of the potential after only a few efolds of infla- 
tion thereby freezing the running of the conformal factor. 
In any case, comparing Eq. (|104f> to the current solar- 
system and astrophysical bounds in Eqs. I|97|) and l|99|) 
shows that scalar-tensor inflation should be considered in 
the search of scalar-tensor effects in cosmology. 



V. CONCLUSION 

In this paper, we have studied the CMB signatures a 
minimal realistic boundary inflation model may have in 
a (reasonably) restricted parameter space in which the 
observable effects come only from the shapes of the pri- 
mordial scalar and tensor power spectra. 

Firstly, it has been shown that the background evo- 
lution involves three different eras corresponding to the 
domination of one of the three fields over the others. Ac- 
cording to the values of the coupling constants, these 
eras may or may not be of inflationary kind. In our 
case, one of the moduli field, namely ip, is associated to 
a non-accelerated expansion during which it rapidly re- 
laxes toward vanishing values. Once ip trapped into its 
minimum, two smoothly connected inflationary periods 
driven by the moduli ip and the matter field \ take place 
until the matter field oscillates around the minimum of 
its potential thereby starting a reheating period. 

At the perturbation level, we have discussed the gen- 
eration and observability of the primordial cosmologi- 
cal perturbations arising during these inflationary eras. 
Since the ^-dominated era precedes the inflationary eras, 
this field has not significant observable effects. Moreover, 
the existence of an attractor in the background fields evo- 
lution during inflation erases any memory of the initial 
conditions when one is concerned with the observable 
perturbations. As a result, only the effective mass m e ff of 
the matter field and the moduli coupling constant c p end 
up being of observable interest for the CMB. The former 
encodes the amplitude of the primordial perturbations at 
Hubble exit and the latter quantifies the changing rate of 
the conformal factor with respect to the evolution of the 
moduli tp. Moreover, as expected for a multifield system, 
there are entropy perturbations which source the adia- 
batic modes after Hubble exit. However, for cosmologi- 
cally relevant values of the moduli coupling constant c p , 
the entropy and adiabatic modes are found to be weakly 
coupled and the adiabatic power spectrum dominates at 
the end of inflation. Nevertheless, even in this regime, the 
running of the conformal factor at Hubble exit has been 
shown to significantly increase the tilt of the adiabatic 
power spectrum. 

We have then assumed a toy cosmological framework 
in order to compute the seeded CMB anisotropies. In this 
framework, the moduli ip and ip survive in the late-time 
cosmology whereas the standard cosmological fluids are 
produced by the decay of the matter field x- I n order to 
study the observable effects stemming from inflation only, 
we have focused on values of the moduli coupling con- 
stant which do not strongly modify the late-time cosmol- 
ogy, namely c p < 0.2. Under this prior choice, we have 
computed the induced CMB anisotropies and performed 
a MCMC analysis of the parameter space given the first 
year WMAP data. The boundary inflation model ap- 
pears to be indistinguishable from a single field chaotic 
model as long as c p < 10 -3 whereas it is disfavored for 
larger values of the coupling constant. Since the overall 
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best fit model lies in the c p < 10 -3 region of the param- 
eter space, the boundary inflation model is statistically 
disfavored compared to a single chaotic model by its two 
additional degrees of freedom which do not help to im- 
prove the fit. Moreover, the current WMAP data lead to 
a 95% marginalized probability bound log(c p ) < —1.64 
which corresponds to the post-Newtonian Eddington pa- 
rameter upper limit 1 — 7 < 2 x 10~ 3 . 

The above bound is not competitive with the solar- 
system upper limit and remains only slightly stronger 
than the late-time cosmological one. However, it holds 
at the time of inflation and provides in this respect a 
very early constraint to the scalar-tensor models which 
would behave during inflation as the boundary model, 
but would relax toward General Relativity afterwards. A 
word of caution is in order since this constraint certainly 
does not apply to all scalar-tensor inflationary models. 
Indeed, it essentially relies on the running of the confor- 
mal factor during the generation of the primordial per- 
turbations, and also depends on the shape of the matter 
potential. For example, one may imagine to freeze the 
running of the conformal factor by stabilising the moduli 
in their bulk potential [see Eq. JSJ]. However, in view of 
the future more accurate data, this work suggests that a 
fully consistent derivation of scalar-tensor theory bounds 
in the context of cosmology should involve both an in- 
flationary and post-inflationary modelisation of the ex- 
pected deviations from General Relativity. 

The present work has been focused on the inflation- 
ary eras only and for simplicity we have not considered 



reheating phenomena that might modify the evolution 
of the cosmological perturbations. It could be interest- 
ing in future works to quantify these effects in view of the 
CMB data. In addition, we have not considered the post- 
inflationary scalar-tensor effects. A natural extension of 
the current work would be to perform a MCMC explo- 
ration on the moduli coupling constants by using simul- 
taneously scalar-tensor inflationary and post-inflationary 
codes to compute the cosmological perturbations. In par- 
ticular, one could quantify the late-time effects of the 
entropy perturbations existing generically between the 
moduli and the cosmological fluids in such brane world 
scenarios. 



Acknowledgments 

We thank Nicole Audiffren, Carlo Contaldi, Gillcs 
Esposito-Farese, Samuel Leach, Jerome Martin, 
Christophe Rhodes, Roberto Trotta and Jean-Philippe 
Uzan for enlightening discussions and advices on the 
different parts of this work. The computations have 
been performed for one part in the Cen tre I nformatique 
National de l'Enseignement Superieur |l22j during one 
CPU-time year, and for the other thanks to a substantial 
time allocation by t he French data processing center 
for Planck-HFI [l23| and by the U. K. Computational 
Cosmology Consortium |l24j | . This work is supported in 
part by PPARC. 



[1] S. Eidelman, K. Hayes, K. Olive, M. Aguilar-Benitez, 
C. Amsler, D. Asner, K. Babu, R. Barnett, J. Beringer, 
P. Burchat , et al., Physics Letters B 592, 1+ (2004), 
URL |http : // pdg . lbl . gov| 

[2] G. Nordstrom, Phys. Zeitschr. 15, 504 (1914). 

[3] T. Kaluza, Sitzungsber. Preuss. Akad. Wiss. Berlin p. 
966 (1921). 

[4] O. Klein, Z. Phys. 37, 895 (1926). 

[5] L. ORaifeartaigh and N. Straumann, Rev. Mod. Phys 
72, 1 (2000). 

[6] J. Polchinski, String theory. An introduction to the 
bosonic string, Vol. I (Cambridge University Press, 
Cambridge, UK, 1998). 

[7] J. Polchinski, String theory. Super string theory and be- 
yond, Vol. II (Cambridge University Press, Cambridge, 
UK, 1998). 

[8] P. Horava and E. Witten, Nucl. Phys. B475, 94 (1996), 

hep-th/9603142. 
[9] E. Witten, Nucl. Phys. B471, 135 (1996), hep- 

th/9602070. 

[10] D. Langlois, Prog. Theor. Phys. Suppl. 148, 181 (2003), 

hep-th/0209261. 
[11] R. Maartens, Living Rev. Rel. 7, 7 (2004), gr- 

qc/0312059. 

[12] P. Brax, C. van de Bruck, and A.-C. Davis, Rept. Prog. 

Phys. 67, 2183 (2004), hep-th/0404011. 
[13] K. Akama, Lect. Notes Phys. 176, 267 (1982), hep- 



th/0001113. 

[14] V. A. Rubakov and M. E. Shaposhnikov, Phys. Lett. 

B125, 136 (1983). 
[15] A.-C. Davis, S. C. Davis, W. B. Perkins, and I. R. Ver- 
non, Phys. Lett. B504, 254 (2001), hep-ph/0008132. 
[16] R. A. Battye and B. Carter, Phys. Lett. B357, 29 

(1995), hep-ph/9508300. 
[17] R. A. Battye, B. Carter, A. Mennim, and J. -P. Uzan, 

Phys. Rev. D64, 124007 (2001), hep-th/0105091. 
[18] P. Binetruy, C. Deffayet, and D. Langlois, Nucl. Phys. 

B565, 269 (2000), hep-th/9905012. 
[19] R. Maartens, D. Wands, B. A. Bassett, and I. Heard, 

Phys. Rev. D62, 041301 (2000), hep-ph/9912464. 
[20] L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 4690 

(1999), hep-th/9906064. 
[21] L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 3370 

(1999) , hep-ph/9905221. 

[22] J. Garriga and T. Tanaka, Phys. Rev. Lett. 84, 2778 

(2000) , hep-th/9911055. 

[23] N. Arkani-Hamed, S. Dimopoulos, and G. R. Dvali, 
Phys. Lett. B429, 263 (1998), hep-ph/9803315. 

[24] I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos, and 
G. R. Dvali, Phys. Lett. B436, 257 (1998), hep- 
ph/9804398. 

[25] C. Ringeval, T. Boehm, and R. Durrer (2003), hep- 
th/0307100. 

[26] G. R. Dvali, G. Gabadadze, and M. Porrati, Phys. Lett. 



23 



B485, 208 (2000), hep-th/0005016. 
[27] I. Antoniadis, R. Minasian, and P. Vanhove, Nucl. Phys. 

B648, 69 (2003), hep-th/0209030. 
[28] D. Cremades, L. E. Ibanez, and F. Marchesano, Nucl. 

Phys. B643, 93 (2002), hep-th/0205074. 
[29] C. Kokorelis, Nucl. Phys. B677, 115 (2004), hep- 

th/0207234. 

[30] E. Kohlprath, Nucl. Phys. B697, 243 (2004), hep- 
th/0311251. 

[31] M. Kolanovic, M. Porrati, and J.-W. Rombouts, Phys. 

Rev. D68, 064018 (2003), hep-th/0304148. 
[32] C. Ringeval and J.-W. Rombouts, Phys. Rev. D71, 

044001 (2005), hep-th/0411282. 
[33] M. Cvetic and H. H. Soleng, Phys. Rept. 282, 159 

(1997), hep-th/9604090. 
[34] F. Bonjour, C. Charmousis, and R. Gregory, Class. 

Quant. Grav. 16, 2427 (1999), gr-qc/9902081. 
[35] N. D. Antunes, E. J. Copeland, M. Hindmarsh, and 

A. Lukas (2002), hep-th/0208219. 
[36] C. Ringeval, P. Peter, and J.-P. Uzan, Phys. Rev. D71, 

104018 (2005), hep-th/0301172. 
[37] K. Koyama, D. Langlois, R. Maartens, and D. Wands, 

JCAP 0411, 002 (2004), hep-th/0408222. 
[38] C. Cartier and R. Durrer, Phys. Rev. D71, 064022 

(2005), hep-th/0409287. 
[39] A. Lukas, B. A. Ovrut, and D. Waldram, Phys. Rev. 

D61, 023506 (2000), hep-th/9902071. 
[40] A. Lukas, B. A. Ovrut, K. S. Stelle, and D. Waldram, 

Nucl. Phys. B552, 246 (1999), hep-th/9806051. 
[41] P. Brax and A. C. Davis, Phys. Lett. B497, 289 (2001), 

hep-th/0011045. 
[42] P. R. Ashcroft, C. van de Bruck, and A. C. Davis, Phys. 

Rev. D69, 083516 (2004), astro-ph/0210597. 
[43] P. Brax, C. van de Bruck, A. C. Davis, and C. S. Rhodes, 

Phys. Rev. D67, 023512 (2003), hep-th/0209158. 
[44] S. Kobayashi and K. Koyama, JHEP 12, 056 (2002), 

hep-th/0210029. 
[45] P. Brax, C. van de Bruck, A. C. Davis, and C. S. Rhodes 

(2003) , hep-ph/0309180. 

[46] P. R. Ashcroft, C. van dc Bruck, and A. C. Davis, Phys. 

Rev. D69, 063519 (2004), astro-ph/0310643. 
[47] S. Kanno and J. Soda, Gen. Rel. Grav. 36, 689 (2004), 

hep-th/0303203. 
[48] G. A. Palma and A.-C. Davis, Phys. Rev. D70, 064021 

(2004) , hep-th/0406091. 

[49] A. C. Davis, P. Brax, and C. van de Bruck, Nucl. Phys. 

Proc. Suppl. 148, 64 (2005), astro-ph/0503467. 
[50] T. Damour and G. Esposito-Farese, Class. Quant. Grav. 

9, 2093 (1992). 

[51] C. M. Will, Living Rev. Rel. 4, 4 (2001), gr-qc/0103036. 
[52] T. Damour and B. Pichon, Phys. Rev. D59, 123502 

(1999), astro-ph/9807176. 
[53] F. Perrotta, C. Baccigalupi, and S. Matarrese, Phys. 

Rev. D61, 023507 (2000), astro-ph/9906066. 
[54] X.-l. Chen and M. Kamionkowski, Phys. Rev. D60, 

104036 (1999), astro-ph/9905368. 
[55] C. Baccigalupi, S. Matarrese, and F. Perrotta, Phys. 

Rev. D62, 123510 (2000), astro-ph/0005543. 
[56] L. Amendola, Phys. Rev. Lett. 86, 196 (2001), astro- 

ph/0006300. 

[57] A. Riazuelo and J.-P. Uzan, Phys. Rev. D66, 023525 

(2002), astro-ph/0107386. 
[58] R. Nagata, T. Chiba, and N. Sugiyama, Phys. Rev. 

D69, 083512 (2004), astro-ph/0311274. 



[59] V. Acquaviva, C. Baccigalupi, S. M. Leach, A. R. Lid- 
die, and F. Perrotta, Phys. Rev. D71, 104025 (2005), 
astro-ph/0412052. 

[60] G. Esposito-Farese (2004), gr-qc/0402007. 

[61] J.-P. Uzan, AIP Conf. Proc. 736, 3 (2005), astro- 
ph/0409424. 

[62] C. Schimd, J.-P. Uzan, and A. Riazuelo, Phys. Rev. 

D71, 083512 (2005), astro-ph/0412120. 
[63] C. S. Rhodes, C. van de Bruck, P. Brax, and A. C. Davis, 

Phys. Rev. D68, 083511 (2003), astro-ph/0306343. 
[64] T. Damour and K. Nordtvedt, Phys. Rev. D48, 3436 

(1993) . 

[65] A. Kogut et al., Astrophys. J. Suppl. 148, 161 (2003), 

astro-ph/0302213. 
[66] L. Verde et al., Astrophys. J. Suppl. 148, 195 (2003), 

astro-ph/0302218. 
[67] G. Hinshaw et al., Astrophys. J. Suppl. 148, 135 (2003), 

astro-ph/0302217. 
[68] G. A. Palma and A.-C. Davis, Phys. Rev. D70, 106003 

(2004) , hep-th/0407036. 

[69] A. Lewis, A. Challinor, and A. Lasenby, Astrophys. J. 

538, 473 (2000), astro-ph/9911177. 
[70] A. Lewis and S. Bridle, Phys. Rev. D66, 103511 (2002), 

astro-ph/0205436. 
[71] L. Kofman et al., JHEP 05, 030 (2004), hep- 

th/0403001. 

[72] N. A. Koshelev, Grav. Cosmol. 10, 289 (2004), astro- 
ph/0501600. 

[73] C. B. Collins, Commun. Math. Phys. 23, 137 (1971). 
[74] C. B. Collins, Commun. Math. Phys. 27, 37 (1972). 
[75] I. S. Shikin, Sov. Phys. JETP 36, 811 (1973). 
[76] C. Gordon, D. Wands, B. A. Bassett, and R. Maartens, 

Phys. Rev. D63, 023506 (2001), astro-ph/0009131. 
[77] F. Di Marco, F. Finelli, and R. Brandenberger, Phys. 

Rev. D67, 063512 (2003), astro-ph/0211276. 
[78] G. Esposito-Farese and D. Polarski, Phys. Rev. D63, 

063504 (2001), gr-qc/0009034. 
[79] A. Linde, JHEP 11, 052 (2001), hep-th/0110195. 
[80] H. Noh and J.-c. Hwang, Phys. Lett. B515, 231 (2001), 

astro-ph/0107069. 
[81] F. Di Marco and F. Finelli, Phys. Rev. D71, 123502 

(2005) , astro-ph/0505198. 

[82] J. Martin (2004), hep-th/0406011. 

[83] M. S. Turner, Phys. Rev. D28, 1243 (1983). 

[84] Y. Shtanov, J. H. Traschen, and R. H. Brandenberger, 

Phys. Rev. D51, 5438 (1995), hep-ph/9407247. 
[85] L. Kofman, A. D. Linde, and A. A. Starobinsky, Phys. 

Rev. D56, 3258 (1997), hep-ph/9704452. 
[86] S. Groot Nibbelink and B. J. W. van Tent, Class. Quant. 

Grav. 19, 613 (2002), hep-ph/0107272. 
[87] J.-c. Hwang and H. Noh, Class. Quant. Grav. 19, 527 

(2002), astro-ph/0103244. 
[88] V. N. Lukash, Sov. Phys. JETP 52, 807 (1980). 
[89] D. H. Lyth, Phys. Rev. D31, 1792 (1985). 
[90] H. Kodama and M. Sasaki, Prog. Theor. Phys. Suppl. 

78, 1 (1984). 

[91] V. F. Mukhanov, H. A. Feldman, and R. H. Branden- 
berger, Phys. Rept. 215, 203 (1992). 
[92] R. Durrer, Fundamentals of Cosmic Physics 15,3, 209 

(1994) , astro-ph/9311041. 

[93] J. Martin and D. J. Schwarz, Phys. Rev. D57, 3302 

(1998), gr-qc/9704049. 
[94] B. Carter (1997), hep-th/9705172. 

[95] D. Langlois, Phys. Rev. D59, 123512 (1999), astro- 



24 



ph/9906080. 

D. S. Salopek, J. R. Bond, and J. M. Bardeen, Phys. [110 
Rev. D40, 1753 (1989). 

S. Tsujikawa, D. Parkinson, and B. A. Bassett, Phys. [Ill 
Rev. D67, 083516 (2003), astro-ph/0210322. 
J. Martin and C. Ringeval, Phys. Rev. D69, 083515 [112 
(2004), astro-ph/0310382. 

R. Easther, W. H. Kinney, and H. Peiris (2005), astro- [113 
ph/0505426. 

N. Kaloper, M. Kleban, A. Lawrence, S. Shenker, and [114 
L. Susskind, JHEP 11, 037 (2002), hep-th/0209231. 
J. Martin and R. Brandenberger, Phys. Rev. D68, [115 
063513 (2003), hep-th/0305161. 

A. R. Liddle and D. H. Lyth, Phys. Rept. 231, 1 (1993), [116 
astro-ph/9303019. 

J. A. Adams, G. G. Ross, and S. Sarkar, Nucl. Phys. 
B503, 405 (1997), hep-ph/9704286. [117 
N. Kaloper and M. Kaplinghat, Phys. Rev. D68, 123522 [118 
(2003), hep-th/0307016. 

P. Hunt and S. Sarkar, Phys. Rev. D70, 103518 (2004), [119 
astro-ph/0408138. 

J. Martin and D. J. Schwarz, Phys. Rev. D62, 103520 [120 
(2000), astro-ph/9911225. 

A. A. Starobinsky and J. Yokoyama (1994), gr- [121 
qc/9502002. [122 
J. Garcia-Bellido and D. Wands, Phys. Rev. D52, 6739 [123 
(1995), gr-qc/9506050. [124 
A. R. Liddle and S. M. Leach, Phys. Rev. D68, 103503 



(2003), astro-ph/0305263. 

B. Bertotti, L. less, and P. Tortora, Nature 425, 374 
(2003). 

J. G. Williams, X. X. Newhall, and J. O. Dickey, Phys. 
Rev. D53, 6730 (1996). 

M. Tegmark and M. Zaldarriaga, Phys. Rev. D66, 
103508 (2002), astro- ph/0207047. 

H. V. Peiris et al., Astrophys. J. Suppl. 148, 213 (2003), 
astro-ph/0302225. 

S. M. Leach and A. R. Liddle, Phys. Rev. D68, 123508 

(2003) , astro-ph/0306305. 

V. Barger, H.-S. Lee, and D. Marfatia, Phys. Lett. 
B565, 33 (2003), hep-ph/0302150. 

J. R. Bond, C. R. Contaldi, A. M. Lewis, and 
D. Pogosyan, Int. J. Theor. Phys. 43, 599 (2004), astro- 
ph/0406195. 

S. Leach (2005), astro-ph/0506390. 

D. N. Spergel et al. (WMAP), Astrophys. J. Suppl. 148, 

175 (2003), astro-ph/0302209. 

M. Tegmark et al. (SDSS), Phys. Rev. D69, 103501 

(2004) , astro-ph/0310723. 

G. Lazarides, R. R. de Austri, and R. Trotta, Phys. Rev. 

D70, 123527 (2004), hep-ph/0409335. 

R. Trotta (2005), astro-ph/0504022. 

URL http://www.cines.fr 

URL http : / / www . planck . f r| 

URL http : //www.damtp . cam. ac .uk/ cosmos 



